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RESUMEN 

El problema restringido de los tres cuerpos es importante en la dinamica de las estrellas dobles y multiples 
y de los sistemas planetarios. Extendemos la version clsica de este problema a una situation que incluye un 
anillo. Encontramos puntos de equilibrio y curvas muy distintas a las del caso clasico. Calculamos el valor del 
exponente de Liapunov para algunas orbitas. 

ABSTRACT 

The restricted three body problem is well-known and very important for dynamics of binary, multiple stars 
and also planetary systems. We extend the classical version of this problem to the situation that there are 
some external forces from the belt. We find that both the equilibrium points and solution curves become quite 
different from the classical case. We also determine the values of Lyapunov Exponent for some important 
orbits. 
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1. INTRODUCTION 

The three body problem is one of the most im- 
portant problems of celestial mechanics and has been 
analytically and numerically studied for centuries. In 
addition to that, three body interaction also plays 
an essential role for dynamics of binary and multiple 
stars. Please see Valtonen (2004) and Dvorak (2004) 
and also their references. 

On the other hand, because there are Asteroid 
Belt and Kuiper Belt for the Solar system, discs of 
dust for extrasolar planetary systems and also cir- 
cumbinary rings for binary systems, these belt-like 
structure should influence the dynamical evolution of 
these systems. For instance, Jiang & Ip (2001) show 
that the origin of orbital elements of the planetary 
system of upsilon Andromedae might be influenced 
by the belt interaction initially. Moreover, Yeh & 
Jiang (2001) studied the orbital migration of scat- 
tered planets. They completely classify the parame- 
ter space and solutions and conclude that the eccen- 
tricity always increases if the planet, which moves on 
circular orbit initially, is scattered to migrate out- 
ward. Thus, the orbital circularization must be im- 
portant for scattered planets if they are now moving 
on nearly circular orbits. 

Therefore, Jiang & Yeh (2003) did some analysis 
on the solutions for dynamical systems of planet-belt 
interaction. In this paper, we further study the effect 
of belts for dynamical evolution of a binary system. 
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2. THE MODEL 

We consider the motion of a test particle influ- 
enced by the gravitational force from the central bi- 
nary and the circumbinary belt. The circumbinary 
belt also provides the frictional force for the test par- 
ticle. 

We assume that two masses of the central binary 
are mi and mi and choose the unit of mass to make 
G(mi + m 2 ) = 1. If we define that 



m 2 



mi + mi ' 



then the two masses are = Gm\ = 1-/2 and 
fi2 = Gm 2 = p.. The separation of central binary is 
set to be unity and fj,± = /j,2 = 0.5 for all numerical 
results in this paper. 

The equation of motion of this problem is (Mur- 
ray & Dermott 1999) 
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where the potential U* is 
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ri = y/(x + fJ.2) 2 + y 2 and r 2 = y/(x- /ii) 2 + y 2 . V 
is the potential from the belt. The belt is a annulus 
with inner radius r,- and outer radius r„, where r, 
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and r a are assumed to be constants. We arbitrary 
set ri — 0.2 and r a — 1.0 for all results in this paper. 

The density profile of the belt is p(r) = c/r p , 
where r = \J x 1 + y 2 , c is a constant completely de- 
termined by the total mass of the belt and p is a 
natural number. In this paper, we set p = 2 for all 
numerical results. Hence, for p = 2, the total mass 
of the belt is 



M h = 



2tt 



p(r')r' dr' d<p = 27rc(ln r a — lnr^). (3) 



The gravitational force /& from the belt is 



Mr) = ~ W 



p(r')r' 
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dr', 
"(4) 

where F(£) and E(£) are elliptic integral of the first 
kind and the second kind. Hence, 
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where /& is in Eq. 0. 

The frictional force should be proportional to the 
surface density of the belt and the velocity of the 
particle. In the x direction, the frictional force is 



-ap(r) 



dx 
~dt 



(G) 



and in the y direction, the frictional force is 

fay = -ap{r)^-, (7) 



where a is the frictional parameter. 

We substitute Eq. © and Eq. 0-0 into Eq.0 
and have the following system: 



% = 2v + x 
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3. EQUILIBRIUM POINTS 

The equilibrium points (x e , y e ) of System 
isfy the following equations 

f(x,y)^x-^ {x + ^ -^-^ 



sat- 





Fig. 1. The curves of f(x, y) = and g(x, y) = (see the 
text for details). 



E 



dr' = 0,(9) 



p(r' 
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dr' = 0. (10) 



In Fig.l, we plot the curves of f(x,y) — (full cir- 
cles) and g(x,y) = (open triangles) for different 
values of Mb. Equilibrium points (x e ,y e ) are inter- 
section of these two. In Fig.l (a), we set Mb = 0, 
so there is no influence from the belt. We find there 
are five usual Lagrangian points, L\, L2, £3, £4 and 
L5 for this case. In Fig. 1 (b), Mb = 0.15 and 
we still have five usual Lagrangian points. In addi- 
tion to that, there are two new equilibrium points 
in the upper half-plane and another two in the lower 
half-plane. In Fig. 1 (c)-(d), we set Mb = 0.3 and 
Mb = 0.5 individually. We also find that in the up- 
per half-plane, there are two new equilibrium points, 
F a and F b . 

4. LYAPUNOV EXPONENT 

Since we have discovered two new equilibrium 
points near L4 (and another two near L5), it would 
be interesting to investigate the orbital behavior 
around these new equilibrium points F a and Fb- For 
a complicated system like ours, it is difficult to rig- 
orously prove if the orbits are chaotic near F a and 
Fb. Nevertheless, we use the calculation of Lyapunov 
Exponent for some orbits whose initial conditions 
are chosen to be close to F a and Fb to understand 
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how sensitively dependent on the initial conditions 
for these orbits. This is in fact one of the most 
important methods to study chaotic systems. We 
follow Wolf et al. (1985) to calculate the values of 
Lyapunov Exponent numerically. To check if our cal- 
culation is correct, we have reproduced the results of 
a given system in their paper. 

In general, the larger value of Lyapunov Expo- 
nent means more sensitively dependent on the ini- 
tial conditions. We choose the initial conditions 
of orbits to be close to the equilibrium points F a , 
F b , L 4 and L 2 individually. Thus, there are 4 dif- 
ferent initial conditions for the orbital calculations. 
To understand the effect of the belt with differ- 
ent mass, we did calculations for 4 different masses 
of the belt (M b = 0,0.15,0.3,0.5) for each chosen 
initial conditions. Although there are no equilib- 
rium points F a , Fb when there is no belt (Mb = 
0), and the locations of equilibrium point F a , Fb, 
L 4 and L 2 would be slightly different for different 
masses of the belt, we still call the initial condition 
(x,y,u,v) = (0.01,0.0225,0,0) Initial Condition F a , 
(x,y,u,v) = (0.01,0.06,0,0) Initial Condition F b , 
(x,y,u,v) = (0.01,1,0,0) Initial Condition L4 and 
(x,y,u,v) — (1.35,0,0,0) Initial Condition L 2 . 

Fig.2(a)-(d) are the results of Lyapunov Expo- 
nent for Initial Condition F a , F b , L4 and L 2 indi- 
vidually. There are 4 curves in each panel of Fig. 2, 
where solid curve is the result of M b = 0, dotted 
curve is the result of M b = 0.15, dashed curve is the 
result of Mb = 0.3 and long dashed curve is the result 
of Mb = 0.5. It is obvious that the values of Lya- 
punov Exponent for Initial Condition F a , Fb is much 
larger than the ones for Initial Condition L4 and L 2 . 
From Panel (a) and (b), we can also see that the 
Lyapunov Exponent for M b — 0.5 and Mb = 0.3 are 
larger than the values for Mb = 0.15 and M b = 0. 
Interestingly, for orbits with Initial Condition L4, 
the values of Lyapunov Exponent are slightly larger 
for Mb = as we can see in Fig. 2(c). In general, 
their values are small for both Initial Condition L4 
and L 2 . The values of Lyapunov Exponent for orbits 
with Initial Condition L 2 approach to when t tends 
to infinite. The orbits are obviously not chaotic for 
this case. 

5. ORBITS 

In this section, we will discuss all the orbits whose 
results of Lyapunov Exponent have been shown and 
discussed in last section. 

Fig. 3, 4, 5 and 6 are the orbits on x — y plane 
for Initial Condition F a , Fb, L4 and L 2 . There are 4 
panels for each figure. Panel (a) is the result when 
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Fig. 2. Lyapunov Exponent (see the text for details). 

there is no belt, i.e. Mb = 0, panel (b) is the result 
when Mb = 0.15, panel (c) is for M b = 0.3 and panel 
(d) is the result for M b = 0.5. 

If one looks at all these 4 figures of orbits at the 
same time, one can immediately understand that it 
seems the orbits with Initial Condition F a and F b 
are much more chaotic than the orbits with Initial 
Condition L4 and L 2 . This impression is completely 
consistent with the one we get from the values of 
Lyapunov Exponent. 

To compare Fig. 3(a) with Fig. 4(a), we found that 
the orbits are similar for Initial Condition F a and F b . 
However, from the comparison between Fig. 3(b) and 
Fig. 4(b), we found that the orbits are quite different 
for Initial Condition F a and F b . These two com- 
parisons show that the existence of belt does make 
the orbits become more sensitive to the initial con- 
ditions. 

6. CONCLUDING REMARKS 

We have provided the equations for a model 
which modifies the classical restricted three body 
problem by including the influence from a belt 
around the central binary. We found that, in addi- 
tion to the usual Lagrangian points, there are two 
new equilibrium points, which we call F a and F b 
around L4 (similarly, there are another two new equi- 
librium points close to L5). 

To study the orbits around these new equilibrium 
points, we calculate the values of Lyapunov Expo- 
nents for orbits with 4 different initial conditions, 
which are close to F a and F b , L4 and L 2 individu- 
ally. We found that the belt makes the system even 
more sensitive to the initial conditions for the orbits 
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Fig. 3. The orbits with Initial Condition F a 



Fig. 5. The orbits with Initial Condition L4. 
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Fig. 4. The orbits with Initial Condition Fb- 

with Initial Condition F a and but does not make 
too much different for the orbits with Initial Condi- 
tion L4 and Z<2- Because the equilibrium points F a 
and Fb happen to be near inner part of the belt and 
Lagrangian points L4 and L2 happen to be around or 
out of outer part of the belt in our system, it seems 
that the orbits near inner part of the belt might be 
more unpredictable than the ones around outer part. 



Fig. 6. The orbits with Initial Condition Li. 
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